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Abstract 



The saturation of QCD chiral sum rules is reanalyzed in view of the new and 
complete analysis of the ALEPH experimental data on the difference between 
vector and axial- vector correlators (V-A). Ordinary finite energy sum rules 
(FESR) exhibit poor saturation up to energies below the tau-lepton mass. A 
remarkable improvement is achieved by introducing pinched, as well as mini- 
mizing polynomial integral kernels. Both methods are used to determine the 
dimension d = 6 and d = 8 vacuum condensates in the Operator Product Ex- 
pansion, with the results: C>6(2.6 GeV^) = -(0.00226 ± 0.00055) GeV^ , and 
08(2.6 GeV^) = -(0.0054 ±0.0033) GeV® from pinched FESR, and compatible 
values from the minimizing polynomial FESR. Some higher dimensional con- 
densates are also determined, although we argue against extending the analysis 
beyond dimension d = 8. The value of the finite remainder of the (V-A) correla- 
tor at zero momentum is also redetermined: n(0) = —4 Liq — 0.02579±0. 00023. 
The stability and precision of the predictions are significantly improved com- 
pared to earlier calculations using the old ALEPH data. Finally, the role and 
limits of applicability of the Operator Product Expansion in this channel are 
clarified. 
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1 Introduction 



More than twenty five years ago, Shifman, Vainshtein and Zakharov ITl proposed 
to use the Operator Product Expansion (OPE) in hadronic current-current cor- 
relators to extend asymptotic predictions of QCD to low energies. In this ap- 
proach there appear universal vacuum expectation values of quark and gluon 
fields, the so-called vacuum condensates, which have to be extracted from ex- 
periment. This extraction is usually carried out by using methods based on 
dispersion relations. Ultimately, one has to relate error afflicted data in the 
time-like region to asymptotic QCD in the space-like region. Unfortunately, 
this task of analytic continuation constitutes, mathematically, a so-called ill- 
posed problem. In fact, extracting condensates from data is highly sensitive 
to data errors. Not surprisingly, results from different collaborations have not 
been always consistent j^J. The main reason for these inconsistencies was fre- 
quently the impossibility of estimating reliably the errors in the method. With 
the release of the final analysis of the precise measurements of the vector (V) 
and axial-vector (A) spectral functions obtained from tau-lepton decay by the 
ALEPH collaboration 3 , an opportunity has been opened to check the validity 
of QCD sum rules and the extraction of condensates in the light-quark sector 
with unprecedented precision. It is therefore appropriate to reanalyze the data, 
taking into account errors and correlations with the least possible theoretical 
bias. In this paper we attempt such a critical and conservative appraisal for 
the interesting case of chiral sum rules. These sum rules involve the difference 
between the vector and the axial- vector correlators ( V-A) , which vanishes iden- 
tically to all orders in perturbative QCD in the chiral limit. In fact, neglecting 
the light quark masses, the (V-A) two-point function vanishes like in the 
space- like region, where the scale 0(300 MeV) is set by the four-quark con- 
densates. The interest in these sum rules is twofold. Apart from describing a 
QCD order parameter, they determine the leading contributions of the matrix 
elements of the electroweak penguin operators 



where Cq is the charge of the quark q. 

In the time-like region, the chiral spectral function pv-A^q^) should vanish for 
sufficiently large = —q^, but judging from the ALEPH data ^ shown in Fig.l 
the asymptotic regime of local duality does not seem to have been reached, i.e. 
the spectral function docs not vanish even for the highest momenta attainable 
in T-decay. 

Under less stringent assumptions, one would hope that at least global duality 
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Figure 1: The ALEPH data on the vector minus axial- vector spectral function 
vs. perturbativc QCD (solid line). 



should hold in the time-like region. In particular, this should be the case for 
the Weinberg- type sum rules which involve the first and second moment 

of the spectral function. Surprisingly, these sum rules also appear to be poorly 
convergent. A likely source of duality violation could be some non-perturbative 
contribution to the correlator (e.g. due to instantons) which falls off exponen- 
tially in the space-like region but oscillates in the time-like region. From Fig.l it 
is obvious that convergence could be improved by incorporating a weight factor 
which would reduce the non-asymptotic contribution to the spectral integral. 
This can be achieved e.g. by considering so called "pinched sum rules" |S] or 
"minimizing polynomial sum rules" '7'. In view of its importance we have cho- 
sen to reanalyze the issue of duality in chiral sum rules on the basis of the new 
ALEPH measurements. Our analysis leads to results showing a significantly 
improved accuracy. 
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2 Finite energy sum rules 



We begin by defining the vector and axial-vector current correlators 

K^i.^) =ijd'x e''- < Q\T{V^{x)V^m\Q > (1) 

^Xil') ^^jd^x e^"^ < 0\TiA,ix)Aim\0 > (2) 

= hgt^i^q^ + qiiiu) n^lg^) - g^g^ no(g^) , 

where V^{x) =: q[x)-i ^_,q{x) :, A^{x) =: g(a;)7^ 759(0;) :, and q (u, d). Here 
we shall concentrate on the chiral correlator Ily.yi = Ily — Ha- This corre- 
lator vanishes identically in the chiral limit (to^ = 0), to all orders in QCD 
perturbation theory. Rcnormalon ambiguities are thus avoided. To define our 
normalization we note that in perturbative QCD 

i imn^^- (.) = i imn^^- (^) - A (1 + - + ••■) (3) 

Non-perturbative contributions due to vacuum condensates contribute to this 
two-point function starting with dimension d = 6, and involving the four-quark 
condensate. The Operator Product Expansion (OPE) of the chiral correlator 
can be written as 



n{Q^)\v-A = E 7W ^?2iv(0',/^') < 02iv(M') > , (4) 



with = —q^. The scale parameter /i separates the long distance non- 
perturbative effects associated with the condensates < 02Ar(/^^) > from the 
short distance effects which are included in the Wilson coefficients C2n{Q^ , f^^)- 
The OPE is valid for complex q^ and moderately large |q^| sufficiently far away 
from the positive real axis. Radiative corrections to the dimension d = 6 contri- 
bution are known They depend on the regularization scheme, implying 
that the value of the condensate itself is a scheme dependent quantity. Explic- 
itly, 
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in the MS scheme, and assuming vacuum saturation of the four-quark conden- 
sate. Radiative corrections for d > 8 are not known. To facilitate comparison 
with current conventions in the hterature it will be convenient to absorb the Wil- 
son coefficients, including radiative corrections, into the operators, and rewrite 
Eq.|01) compactly as 



A'>3 ^ 

where we are dropping the subscript (V-A) for simplicity. We will be concerned 
with Finite Energy Sum Rules (FESR), which are nothing but the Cauchy 
integral 

1 1 f 

(7) 

where PAr(s) is an arbitrary polynomial, i.e. 

Pwis) = ao + ais + a2S^ + . . . + aNS^ , (8) 

/tt = 92.4 ± 0.26 MeV 12 , and v{s) (a(s)) is the vector (axial-vector) spectral 
function measured by ALEPH in tau-decay [3|, normalized to the asymptotic 
value 

v{s)qcd - a{s)QCD = i (l + — + ...) . (9) 



2 

The axial-vector spectral function a(s) does not include the pion pole contribu- 
tion which is added separately. For most purposes one can work in the chiral 
limit m^r = 0, i.e. P/v('7iJ) in Eq.© may be replaced by ap/^. The standard 
FESR follow from the theorem of residues and assuming the Wilson coefficients 
are just numbers, 

(^)^'^+^) O27v(so) = ^ p ds Hs) - ais)] - fl 5ni {N = 1, 2, 3...) , 

(10) 

where the index N has been rearranged for convenience. Strictly speaking, 
Ea. H10|l only holds for the constant terms of the Wilson coefficients. Otherwise 
condensates of lower or higher dimension get mixed when taking into account 
radiative corrections due to the logarithmic terms (ln(/Lt^/(5^) or higher). How- 
ever this mixing of operators of different dimensions occurs only at order in 
a given FESR For dimensional reasons the contribution of the operators 
of higher dimension in Ea. (|10|) vanishes for large Sq, while that of operators 
of lower dimension increases with sq. The latter contribution is particularly 
disturbing for operators of high dimension. As the logarithmic terms of the 
relevant Wilson coefficients are not known (except the one for Og) we will ne- 
glect the contribution of operators of dimension unequal to 2 in Eq. (|10l) . This 
approximation is inherent to every sum rule analysis of the r-data and can only 
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be justified a posteriori by demonstrating that the right hand side of Ea. lfTU|l 
is (almost) constant. We will, however, examine the order of magnitude of the 
mixing to be expected by using Eq.© to estimate the effect of the radiative 
corrections of Oq. 

For = 1, 2 Ea. HlU|l leads to the first two (Finite Energy) Weinberg sum rules, 
while for = 3, 4 the sum rules project out the d = Q and d — % vacuum 
condensates, respectively (notice that in the chiral limit ©2 = = 0). In 
order to check the convergence of the sum rules we consider the first Weinberg 
sum rule 

^^^'""^ ^ i^r ^""^'^ ' "^'^^ " ^^^^ 

Strictly speaking sq — > cxd, but precocious scaling would imply that the sum 
rule should be saturated at moderate values of sq. From Fig. 2, which shows 
Wi(so), one can see that this is clearly not the case, even at the highest energies 
accessible in r-decay. This lack of precocious scaling in the sum rule can be 
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Figure 2: The first Weinberg sum rule, Ea. (|ll|l as a function of sq together with 
/2 = (0.00854 ±0.00005) GeV^ (solid line). 

simply explained by looking at the measured spectral function in Fig.l. If the 
spectral function had reached its approximate asymptotic value (i. e. zero) 
starting, let us say, from s = 2 GeV^ then the spectral integral of Eq. would 
have yielded for all sq >2 GeV^ . This observation shows us a way out of the 
dilemma by turning to the more general sum rules of Eq.Q. One can choose 
the polynomial PAr(s) in the sum rule Eq.Q in such a way that the problematic 
contribution of the integration region near the endpoint of the physical cut is 
minimized. This method addresses two problems at the same time, the first 
being that experimental errors of the spectral functions grow considerably near 
the limit of phase space, and the second, that the asymptotic QCD formula 
is unreliable on the contour region near the physical cut. In following this 
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method we will employ two types of sum rules, pinched FESR and minimizing 
polynomial sum rules. 



3 Pinched sum rules 

We begin by considering a linear combination of the first two Weinberg sum 
rules 

ds{l--)[v{s)-a{s)]^fl. (12) 

The left hand side of Ea. p2|l as a function of so is shown in Fig. 3, together 
with the right hand side. It is very reassuring that the sum rule appears to be 
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Figure 3: Pinched Weinberg sum rule, Ea. (|12|l . as a function of sq. The solid 
line is fl = (0.00854 ± 0.00005) GeV^ 

saturated for sq > 2.3 GeV^. We note that the error band is about a factor 
three smaller than that found in a similar analysis 14 using the old ALEPH or 
OPAL data The influence of the logarithmic dependence of Og in this 

sum rule is about 4 x 10^® GcV^ in the region of the saturation with sq. 

Motivated by this success, we impose the Weinberg sum rules as constraints 
in other pinched finite energy sum rules involving different moments. To be 
precise, we assume that there are no operators of dimension d — 2 nor d — A, 
which is true in the chiral limit, together with the condition that the FESR 
involves factors of (1 — ^) so as to minimize the contribution near the cut. In 
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this way, we write the Das-Mathur-Okubo sum rule in the form 

n(0) = /" - (1 - -r Hs) ~ a{s)\ + ^ , (13) 
47r^ Jo s So So 

where n(0) is the finite remainder of the chiral correlator at zero momentum. It 
is related to Zio, the counter term of the O(p^) Lagrangian of chiral perturbation 
theory ^7], which has been calculated independently, 



n(0)--4Lio = 



> -Fa 



= 0.026 ±0.001 



(14) 



where < > is the electromagnetic mean squared radius of the pion, < > = 
0.439 ±0.008 fm^ 18 , and Fa is the axial- vector coupling measured in radiative 
pion decay. Fa — 0.0058 ± 0.0008 From Fig. 4 we see that this sum rule is 
even more remarkably satisfied. This can be understood by noting that the sum 
rule emphasizes less the high-s region where duality violation competes against 
stability. Numerically, we find from the DM0 sum rule 




Figure 4: The finite remainder of the chiral correlator at zero momentum, 11(0), 
from Ea. H13(l as a function of sq. The solid line is the central value 11 (0) = 
0.02579. 



n(0) = 0.02579 ± 0.00023 , (15) 

a result showing a remarkable accuracy for a strong interaction parameter. In 
this particular sum rule the contribution from the logarithmic term of Oq van- 
ishes. 

Next, we turn our attention to the extraction of the condensates with the help 
of the pinched sum rules. The philosophy of our calculation is threefold, viz. 
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(i) to assume that dimension d = 2 and d = A operators are absent in the OPE 
of the chiral current , (ii) to require that the polynomial projects out only one 
operator of the OPE at a time, and (iii) to require that the polynomial and 
its first derivative vanish on the integration contour of radius |s| = sq. For 
the caveats on point (ii) of this approach see the text above. In this way one 
obtains for iV > 3 the sum rules (ignoring the energy dependence in the Wilson 
coefficients) 



02Niso) = (-1) 



N-1 




— ds[iN-2y 



N-1 



[vis)-a{.s)]^{N~2).s^-'fA iN>3) 



(16) 



Note that there is always a pinch factor (s — sq)'^ in the polynomial. We use once 
again the new ALEPH spectral function and error correlations in this sum rule. 
The crucial point of the extraction of the condensates is a careful inspection 
of the stability of the result with respect to the variation of all parameters in 
the analysis. In our case there is only one parameter, namely the radius Sq. 
This fact contrasts positively with other approaches based on Laplace sum rules 
which involve at least two parameters, and in addition do not project out just 
one single operator, even for sq oo. As for stability, if a meaningful value of 
O2N is to be extracted we would expect the r.h.s. of Eg. 1)16(1 to be constant for 
all So larger than some some critical value . We can call this requirement strong 
stability. It is best discussed on the basis of the figures below which show the 
predictions for various condensates. Figure 5 shows the result for the dimension 
d = 6 condensate. There is an obvious stability region: 2.3 < So(GeV^) < 3 




s (GeV') 



Figure 5: The dimension d = 6 condensate, Oe, from Ea. (|l()|l as a function of 
so- The sohd line is the central value Oe = -0.00226 GeV^. 
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Figure 6: The dimension d = 8 condensate, Os, from Ea. (|18l ) as a function of 
So- The sohd line is the central value Oe = -0.0053 GeV^ 



from where we find 

06(2.7 GeV^) = -(0.00226 ± 0.00055) GeV^ . (17) 

This value is consistent with the one found from the vacuum saturation approx- 
imation O^s = -0.0020 GeV^ from Eq.© with < qq > (sq) = -0.019 GeV^ 
and a{so)/'K = 0.1 , but it is significantly lower than the one found in some ear- 
lier analyses based on the old, incomplete ALEPH data; e.g. Oq — — (0.004 ± 
0.001) GeV^ obtained in 14 using a similar stability criterion as here. The 
contribution of the logarithmic term from the Og coefficient in H17|l . in the re- 
gion of So considered here, is about 8 x 10^® GeV^ and hence neghgible within 
the errors. 

To facilitate a comparison with an alternative type of sum rules discussed in the 
next section, we give the sum rule for Og explicitly 

^8(^o) = " [2sl - 3sls + s^Ms) - a{s)] + 2sl}l . (18) 

(Notice that 2s|] — 3sqS + = (sq - s)^ (s + 2sq)). The result for Og is shown 
in Fig. 6. In spite of the larger errors there is still a distinct region of duality in 
the interval: 2.3 < so(GeV^) < 3 , which yields 

08(2.6 GeV^) = -(0.0054 ± 0.0033) GeV*^ (19) 

Both the sign and the numerical value of this condensate are controversial (see 
e.g. Table 1 in 0). Our result agrees within errors with e.g. that of l^-fTU). 
and that of ^Hl but disagrees in sign with [201, [12], and with the result of 
the minimal hadronic approximation of large Nc [21 • In |20| the sum rules were 
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evaluated at very low values of so, mainly because the data at that time was 
considered to be too inaccurate at higher sq. Fortunately, this state of affairs 
has now changed with the new ALEPH analysis. Eq.© can be used to estimate 
the effect on Og due to the mixing of Oq arising from the logarithmic term. We 
find a correction of about— 3 x 10~* GeV® which is negligible compared to the 
data error in Eq. ifT^ . The results of the sum rules for do , and ©12 as given by 
Ea. H16|l are shown in Figs. 7 and 8, respectively. The strong stability obtained 
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Figure 7: The dimension d = 10 condensate, Oio, from Eq.(^|) as a function 
of Sq. The solid line is the central value Oio = 0.0355 GeV^*'. 

so far is now no longer obvious, and we find at sq ~ 2.5 GeV^ 

Oio(2.5 GeV^) = (0.036 ± 0.014) GeV^° (20) 
012(2.5 GcV^) = -(0.12 ± 0.05) GeV^^ . (21) 

These results, though, should be taken cum grano salis, to wit. Because the 
OPE is an asymptotic series, the upper limit of the integration range, sqi should 
increase as the dimension of the operators increases. The comparison of the 
numerical results for Og, 0%, Oio, and O12 indicates that at a scale of about 
|s| = 1 GeV^ the OPE starts to diverge at dimension d = 10. In addition, the 
problem of the mixing of operators of different dimensions becomes more severe 
for higher dimensional operators. For these reasons we believe that it is rather 
meaningless to extract quantitative results for condensates of dimension higher 
than d = 8 from r-decay spectral functions. 
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Figure 8: The dimension d = 12 condensate, Ou, from Eg. 1)161 ) as a function 
of So- The soHd hue is the central value O12 = —0.12 GeV^^. 



4 Minimizing Polynomial Sum Rules 



The starting point in this analysis is the general sum rule Eq.Q. The poly- 
nomial can be chosen in such a way that the problematic contribution of the 
integration region near the endpoint of the physical cut is minimized. With the 
normalization condition 

Pn (s = 0) = 1, (22) 

we require that the polynomial P/v(s) should minimize the contribution of the 
continuum in the range [c, so] in a least square sense, i.e. 

/So 
s''PN{s)ds^O {k = Q,...N -1) , (23) 

The parameter c can be chosen freely in the interval < c < sq. On the basis 
of the spectral function of Fig.l, a reasonable choice would be 2 GeV^ < c < 
So ~ 3 GeV^. In a sense, these polynomials are a generalization of pinched 
moments, and the P/v(s) will approach (s — so)^ when c — > so. While pinched 
moments eliminate the contribution on the physical real axis at a single point, 
our polynomials tend to eliminate a whole region (from c to so). The degree N 
of the polynomial can be chosen appropriately to project out certain terms in 
the OPE, Eq.©. The polynomials obtained in this way are closely related to 
the Legendre polynomials as follows. Let us introduce the variable 

:^(.)-^^f^-^+.(0), (24) 
(so - c) (so - c) 
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and define the polynomials as 



where Ln{x) are the Legendre polynomials 

We give here only the first few minimizing polynomials 

P, (.) = 3(2.-(.,+^.))'-(.,-cf ^^gj 

3 (so + c) - (so - c) 
^ 5 (2s - (so + c)f - 3 (so - c)" (2s - (sp + c)) ^^g) 
-5 (so + c)^ + 3 (so - cf (so + c) 



If the polynomials are expressed as in Eq. (8), then from Eqs.Q and pU|) there 
follows the sum rule 

aodo + aidi + ... + aNd-N - fl ^ a20e - asOs + ... + (-1)^ ajvC'2Ar+2 , (30) 
where the constants 

dN = -rT dss^ {v{s) - a{s)) (31) 
4vr^ Jo 

are to be determined from the ALEPH data. 

We begin with the Og condensate and obtain from Ea. H3()(l the sum rule 

do + aidi + 02^2 - fl = a20e . (32) 
After substituting P2 from Ea. (|28|l this sum rule becomes 

Oe{so) = ^{sl + 4soc + c^) {d^ - fl) - (so + c) di + • (33) 

In the sequel we choose the initial value sq = 3 GeV^, but will subsequently 
change it in the range 2.5 < so(GeV^) < 3 in order to verify the criterion of 
strong stability. The condensate Oq from Ea. (|33|l is plotted in Fig. 9 as a 
function of c. One can appreciate a stable point for c — 2.5 GeV^. Fixing c at 
this point of minimal sensitivity has been discussed previously in other FESR 
applications, e.g. in [7|. For c — 2.5 GeV^, and so = 3 GeV^ we obtain from 
Eq.lI22|l the result: 

06(3 GeV*') = -(0.0023 ± 0.0013) GeV^ , (34) 
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Figure 9: The condensate Oq as a function of the parameter c. 



which compares well within errors with the previous result from the pinched 
sum rule, Ea. (|17|l . We have tested positively the stability around this point by 
choosing different values of sq in the range 2.5 < So(GeV^) < 3. For instance, 
using So = 2.5 GeV^ we obtain 06(2.5 GeV*') = -(0.00224 ± 0.00046) GeV^ 
Other values of sq in the above range lead to similar results, well in agreement 
within errors with Ea. H34|l . thus satisfying the criterion of strong stability. 

Next, we consider the Og condensate, and use Ea. H30|) to obtain the sum rule 

do + aidi + 02^2 + asds - fl ^ a20e - a^Os ■ (35) 

The presence of Oq in the sum rule for Og can be dealt with in two ways. One 
could insert the numerical value of Oq, e.g. Ea. l|^ . as obtained from its own 
sum rule, or rather substitute the analytic expression of the sum rule itself. 
The latter procedure yields the best possible results in terms of stability and 
accuracy, and leads to the sum rule 

Osiso) = -I [so(so + 2c)2 + c3] + l [3(s2 + c2) + 4soc] di-dg • (36) 

5 10 

Notice the welcome absence of the term involving the second moment, i.e. d2] it 
cancels out when substituting in Ea. (|35|l the sum rule for Oq, Ea. H33() . Choosing 
again the initial value sq = 3 GeV^, one obtains for Og the results shown in 
Fig. 10. One can see again a stability region near c = 2.5 GeV^, leading to the 
result 

08(3 GeV^) = -(0.0048 ± 0.0039) GeV^ (37) 

It is worth mentioning that the polynomial coefficients entering the sum rule 
Eq. H36|l differ significantly from the ones in the corresponding pinched sum rule 
Ea. H18|l . It is therefore reassuring that both results for Og are compatible. 
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To check for strong stability we have, once again, varied sq in the range 
2.5 < so(GeV^) < 3. For sq = 2.5 GeV^ we find Os{2.5 GeV®) = -(0.0056 ± 
0.0024) GeV^, and similar results for other values of sq in the above range. 
Thus, the criterion of strong stability is again satisfied, albeit within large 
errors. Proceeding beyond dimension d = 8 is marred by the same problems 
mentioned at the end of Section 3; the minimizing polynomial FESR do not 
avoid the divergence of the OPE, nor the increasing importance of operator 
mixing. 
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Figure 10: Og as a function of the parameter c. 



5 Conclusion 



The final ALEPH data for the chiral spectral function v(s) — a{s) shows clearly 
that this spectral function has not yet reached its asymptotic form dictated by 
perturbative QCD, i.e. it does not vanish, even at the highest energies attain- 
able in T-decay. If the asymptotic regime had been reached precociously, let us 
say at ~ 2 GeV^, then it would have been straightforward to calculate the 
non-perturbative condensates with the help of the Cauchy Integral. Since this 
is not the case, some method to improve convergence must be applied. We have 
shown that in the framework of FESR this can be done by suitably reducing 
the impact of the high energy region in the dispersive; intc;gral, cither by using 
pinched sum rules or by using minimizing polynomial sum rules. We first used 
the data in a pinched linear combination of the first two Weinberg sum rules 
which follow from the fact that there are no dimension d = 2 and d = 4 operators 
contributing to the chiral correlator to demonstrate the precocious saturation 
of the sum rule and the remarkable effectiveness of the method. Motivated by 
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this success, we determined a number of QCD condensates by making maximal 
use of the fact that there are no dimension d = 2 and d — A operators and 
requiring strong stability for both methods, i.e. we varied the radius sq in the 
Cauchy integral beginning at the end of r-decay phase space and required that 
the condensates calculated from the data should be reasonably constant for all 
So in some finite region including the end of phase space. We do not assume (as 
is done in most FESR calculations) that the dispersive integral vanishes from 
So — > oo. By showing that there is "strong stability" i.e. precocious saturation 
of the FESR we prove that this region contributes only negligibly. It would 
indeed be surprising if the observed stability would disappear for so larger than 
the end of phase space. We do, however, have to make the assumption, inherent 
in all sum rule analyses of r-decay, that unknown O(a^) effects of mixing of 
operators of different dimensions are negligible for the relevant duality range 
2.5 GeV^ ^ So ^ 3 GcV^. We have checked explicitely that in all the cases 
considered in this work this is the situation when one takes into account the 
logarithmic term of the dimension six Wilson coefficient. The results for Oe 
and Os satisfy this strong stability criterion as is best seen from the figures. 
Extraction of higher condensates of dimension d > 10 leave room for interpre- 
tation, but the conclusion that they grow rapidly with dimensionality is rather 
obvious. Together with the increasing importance of operator mixing, it makes 
the extraction of these condensates a difficult exercise. Our result that Oq and 
Og, have the same sign is in conflict with some of the earlier determinations 
based on the incomplete ALEPH data but agrees with others (see e.g. for a 
comparative summary). 
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